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Using single-crystal transistors, we have performed a systematic experimental study of electronic 
transport through oxidized copper/rubrene interfaces as a function of temperature and bias. We 
find that the measurements can be reproduced quantitatively in terms of the thermionic emission 
theory for Schottky diodes, if the effect of the bias-induced barrier lowering is included. Our analysis 
emphasizes the role of the coupling between metal and molecules, which in our devices is weak due 
to the presence of an oxide layer at the surface of the copper electrodes. 



The electronic transport properties of contacts be- 
tween metals and semiconductors are notoriously sensi- 
tive to the microscopic details of the interface between 
the two materials. For inorganic semiconductors such as 
Silicon, the level of material control is sufficient for the 
quantitative analysis of transport in terms of established 
theoretical models based on the concept of Schottky bar- 
rier [![. For organic semiconductors, on the contrary, the 
control of the interfacial properties (e.g., disorder present 
in the molecular material, the electronic coupling be- 
tween metal and molecules, the density of surface states, 
etc.) is poor and no systematic quantitative compari- 
son between experiments and theory has been possible 
so far [2]. As a consequence, our understanding of trans- 
port across metal/organic interfaces is limited and it is 
not even known whether the conventional Schottky the- 
ory developed for inorganic semiconductors works also 
for organic materials. 

Recently, we have demonstrated that the electrical char- 
acteristics of oxidized Cu, Ni, and Co contacts in or- 
ganic single-crystal field-effect transistors (FETs) exhibit 
an excellent level of reproducibility j3(. Here, we ex- 
ploit this reproducibility to perform a systematic study 
of bias and temperature dependent transport across ox- 
idized copper/rubrene interfaces. We find that the elec- 
trical characteristics of these contacts are well described 
quantitatively by the conventional Schottky theory with 
meaningful values for the microscopic parameters, if the 
bias-induced lowering of the Schottky barrier is taken 
into account. We attribute the good agreement -as well 
as the experimental reproducibility- to the presence of 
the oxide layer at the surface of copper, which causes the 
coupling between metal and molecules to be weak. 
The rubrene transistors consist of « lfim thick rubrene 
crystals laminated onto doped Silicon substrates (act- 
ing as a gate) covered by a 200 nm SiC-2 layer, with 
lithographically defined, electron-beam evaporated Cop- 
per electrodes (see Ref.jj] for details). Prior to lamina- 
tion, the substrate is cleaned in an Oxygen plasma to re- 
move possible resist residues from the dielectric surface. 
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FIG. 1: I — V curves measured at room temperature on a 
short channel FET for three different gate voltages (Vg = 
— 10, —20, — 30V), showing the characteristic behavior of con- 
tact dominated devices (i.e., gate voltage independence and 
steep non-linearity at low bias). The top left inset is an optical 
microscope image of one of our devices (the bar corresponds 
to 100/im). The bottom right inset shows the differential 
conductance normalized at the zero-bias value for 26 different 
samples, illustrating the reproducibility of the measurements. 



This also influences the Cu contacts by enhancing the 
oxidation of their surface and possibly changing the ma- 
terial work function. We have therefore performed photo- 
emission spectroscopy measurements on identically oxi- 
dized copper films, and found a value of 5.5 eV (the mea- 
surements were performed in air with the Riken Keiki 
AC-2 system, on films prepared under identical condi- 
tions used for the device assembly, for more details see 
[5j). This suggests that the Fermi level in oxidized cop- 
per should align well with the highest occupied molecular 
orbital of rubrene, also located roughly 5.5 eV below the 
vacuum level 

Transport measurements were performed on sufficiently 



2 



short channel devices, in which the channel resistance 
is negligible with respect to that of the contacts Q Q- 
In these contact-dominated devices, the I-V characteris- 
tics are gate voltage independent and show a steep non- 
linear increase of the source drain-current at low bias 
(Fig 1). To understand the mechanisms of charge injec- 
tion at the interface we studied the temperature and bias 
dependence of the source-drain current, see (Fig. 2). The 
current decreases with decreasing temperature and if we 
plot ln(UT)/T 2 ) vs 1/T - as normally done for Schottky 
diodes [l| - a linear dependence is observed, indicative of 
thermally activated transport (Fig. 3b). The slope corre- 
sponds to the activation energy, which is plotted in Fig. 3a 
as a function of bias. Measurements performed on many 
different devices exhibit similar behavior and the inset of 
Fig. 3a show that the spread in the measured values of 
activation energy is rather small. 

By modeling the devices as two oppositely biased Schot- 




and temperature. To this end we use the conventional 
theory for Schottky diodes modified to include the 
effect of a tunnel barrier at the metal/semiconductor in- 
terface [IJ, which in our devices is due to the presence 
of the CuOa, layer at the metal surface. Mathematically, 
the presence of a tunnel barrier at the interface does not 
modify the functional dependence of the current on bias, 
and can be included as a prefactor (proportional to the 
transmission coefficient of the barrier) in the expression 
of Io. Physically, however, the presence of a barrier is im- 
portant, because it causes a weak coupling between the 
metal and molecules that, as we discuss below, results 
in a small surface density of states at the surface of the 
organic semiconductor. 

Under reverse bias conditions, Jo depends on V because 
of the voltage-induced Schottky barrier lowering. This 
includes the image-charge barrier lowering (the so-called 
Schottky effect; first bracket in Eq.2) and the field- 
induced change in the dipole due to electrons occupy- 
ing surface states 12| (second bracket in Eq.2). The 



expression for the reverse current then reads [l[ 
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a, defined in Eq. (4) below, is a parameter that quantifies 
the effect of surface states (and q is the electron charge). 
The expression for the forward current reads: 



FIG. 2: I — V characteristic of a contact-dominate device, 
measured at Va = — 30U, for different temperatures. The 
inset shows the temperature dependence of the differential 
conductance (obtained by numerical differentiation of the I — 
V curves) normalized to the value at V=0 and T=285K. In 
both graphs, the open symbols represent the experimental 
measurements and the continuous lines the theoretical values 
calculated using two Schottky diodes in series (shown in the 
scheme on the top right), as described in the text. 

tky diodes connected in series (see scheme in Fig. 2), the 
measured I — V curves can be directly related to the 
electrical characteristics of the source and drain contacts. 
Using this approach, we have previously shown that the 
usual expression I(V, T) = Iq ( e eV / nkT — i) for the current 
through a Schottky diode (n ~ 1 is the ideality factor), 
with Iq taken as constant, reproduces the low bias peak 
in the differential conductance (see inset of Fig. 2), but 
not the behavior of the I — V curves at bias higher than 
few times kT [3]. Here, we analyze quantitatively the 
data by taking into account full dependence of Iq on bias 
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Eq. 1-3 are valid in the thermionic emission regime, 
which is expected to be correct for low-transparency con- 
tacts, as discussed in Ref. §. Using these expressions, 
we calculate numerically the I — V curves through the 
two diodes connected in series. 

In order to compare the results of the calculations with 
the experiments, the values of the parameters in Eq.l- 
3 need to be fixed. Specifically, Nry is the density of 
dopants (unintentionally present) in the rubrene crystals, 
recently measured [13[ to be N D w 1 — 5 1 20 m~ 3 . Vbi 
is the built in potential, which is theoretically expected 
to be smaller than the barrier height The precise 
value in our devices is not known and in the calculation 
we fix Vbi = 80 meV (smaller than the measured activa- 
tion energy) . We have checked that small changes in this 
value do not affect the quality of the comparison between 
theory and experiments, es = 3 is the dielectric constant 
of rubrene. / is a constant that determines the absolute 
magnitude of the current, which is proportional to the 
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junction area and to the transmission probability across 
the CuO x layer. It drops out if we confine our analysis 
to the normalized differential conductance (see inset of 
Fig. 2). The height of the Schottky barrier <j> and a are 
used as fitting parameters. Whereas the value of a is not 
known a priori, the value of <E>o can be estimated by ex- 
trapolating at zero bias the activation energy measured 
at high bias (see Fig. 3) and cannot be varied much with- 
out affecting the comparison between theory and data. 
The continuous lines in Fig. 2 and 3 show the results of 
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FIG. 3: In (a) the measured activation energy (open dots) 
is compared with the calculations (continuos line). The in- 
set shows the activation energy measured on several different 
devices on (SiC>2 and Ta2 0s), illustrating the relatively small 
sample-to-sample variation. In (b) the linearity of the plot of 
ln(I/T) vs 1/T indicates that the thermally activated behav- 
ior is experimentally well obeyed for different values of bias 
voltage (the lines are guides to the eye). 

the calculations. Both for the current and for the differ- 
ential conductance in Fig. 2 an excellent agreement with 
the experiments is found. The bias dependence of the 
activation energy is also reproduced satisfactorily, with 
only a small deviation at low bias. This quantitative 
agreement is remarkable, especially if one considers that 
the theoretical expressions that we have used are strictly 
speaking valid for a one-dimensional geometry, whereas 
our transistors are planar devices. The best fits are ob- 
tained for <i>o is 0.13 — 0.15 eV and a ~ 2 nm, with the 
precise values depending on the specific device. 
In order to interpret the value of a we recall that accord- 
ing to theory [l4| 
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Here, D$ is the density of states at the metal/organic in- 
terface (inside the HOMO-LUMO gap) . 5 is the thickness 
of the interfacial dipole layer, which is determined by the 
decay length of the wave-function of these states inside 
the bulk of the semiconductor 14 1 . Its value can be 
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of the Schottky barrier). We obtain 
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where m is the mass of the charge carrier in the HOMO 
band of rubrene, recently found to be close to the free 



electron mass |15(). Taking this value and the value of 
a obtained from fitting the data, we use Eq. (4) to es- 
timate the surface density of states. We find D$ ~ 5 
10 11 - 1 10 12 eV^cm -2 , between one and two orders 
of magnitude lower than what is tipically measured in 
clean metal/organic semiconductor interfaces by means 
of photo-emission spectroscopy 

The small density of surface states is consistent with the 
presence of the CuO^ layer at the surface of the cop- 
per electrodes. In fact, in an organic semiconductor like 
rubrene, the mid-gap surface states are induced by hy- 
bridization of molecular states due to the coupling to 
the nearby metal [l6| (since rubrene molecules are van 
der Waals bonded, no dangling bonds are present at the 
crystal surface). If the rubrene molecules are only weakly 
coupled to the metal -owing to the presence of the oxide 
layer acting as a tunnel barrier- hybridization is strongly 
suppressed and so is the density of surface states 
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Our analysis, therefore, points to the crucial role of the 
coupling between metal and molecules in determining the 
I - V characteristic of metal/organic contacts. The cou- 
pling has two main effects on transport. On the one 
hand, a weak coupling tends to increase the contact re- 
sistance, because charge carriers have to tunnel through 
a low transmission barrier. On the other hand, a weak 
coupling suppresses the formation of surface states and 
prevents the formation of large dipoles, which can sub- 
stantially increase the height of the Schottky barrier. The 
competition between these two mechanisms determines 
the value of the contact resistance (see also [18| for beau- 
tiful recent experiments qualitatevely consistent with this 
conclusion) . 

In summary, we have used organic single-crystal transis- 
tors to show that the conventional theory for transport 
through Schottky barriers does describe quantitatively 
the properties of metal/organic interfaces, at least in the 
case in which the coupling between metal and molecules 
is weak. 
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estimated by considering the exponential damping of the 
wavefunction of an electron under a potential barrier U 
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